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^  ABSTRACT 

Preliminary  efforts  are  made  to  develop  a  methodology  for  analyzing  the 
progressive  failure  of  projectile/target  systems  under  hypervelocity  impact. 
Particular  emphasis  is  placed  on  the  thermal/mechanical  interaction  effects 
not  only  for  evaluating  the  change  in  the  local  strain  rates  and  material 
properties  but  also  in  finding  the  temperature  and  latent  heat  at  which  the 
solid  may  locally  transform  into  liquid  and/or  gas.  Such  a  transition  can 
be  accurately  determined  by  the  inflection  points  on  the  time  dependent 
H-function  curve  that  serves  as  a  measure  of  the  order  and/,  r  disorder  of  a 
physical  system. 

Damage  of  a  tungsten  projectile  with  a  1:3.57  aspect  ratio  impacting  a 
5083  aluminum  target;  with  a  1:3  aspect  ratio  at  9,000  m/sec  is  evaluated. 
Melting  of  the  target  material  is  predicted  in  the  region  directly  underneath 
the  projectile  corner.  Geometrically  and  metal lurgically  similar  systems 
scaled  down  one-half  and  one-quarter  in  size  are  also  considered  resulting 
in  damage  zones  that  are  not  proportional.  Such  information  is  useful  for 
developing  scaling  relations  even  though  the  results  are  nonlinear  functions 
of  the  governing,  parameters . 
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FORWARD 


This  investigation  on  hypervelocity  impact  has  been  conducted  by  the 
Comp tek  Corporation  for  the  Army  Materials  Technology  Laboratories  in 
Watertown,  Massachussetts  under  Contract  DAAG46-85-C-0061 . 

The  principal  investigator  for  the  study  was  Dr.  George  C.  Sih  who 
has  successfully  developed  the  basic  thermodynamics  for  analyzing  material 
transformation  under  high  speed  impact.  Analyzed  in  this  report  is  the 
melting  of  5083  aluminum  target  by  a  tungsten  projectile.  Initial  impact 
results  are  also  obtained  for  projectile/target  systems  in  full  scale, 
half  scale  and  quarter  scale  with  the  objective  of  developing  sub-scaling 
relations. 

Dr.  S.  C.  Chou  has  been  the  AMTL  technical  monitor.  The  author 
wishes  to  thank  him  for  providing  some  of  the  pertinent  information.  The 
encouragement  and  interest  expressed  by  Mr.  John  F.  Dignam  during  the 
course  of  this  project  are  also  gratefully  acknowledged. 
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I.  INTRODUCTION 

When  a  solid  comes  in  contact  with  another  at  high  speed,  there  prevails 
locally  an  extremely  high  rate  of  energy  transfer.  The  initial  kinetic 
energy  of  the  moving  body  is  converted  to  different  forms  depending  on  the 
impact  speed,  body  geometry  and  size  and  material  property.  If  the  impact 
velocity  is  sufficiently  high,  local  material  transformation  can  take  place 
such  that  the  solid  may  melt  or  even  vaporize  in  addition  to  damage  by  per¬ 
manent  deformation  and/or  fracture.  Material  transformation  and  failure 
behavior  are  not  only  time  varying  but  they  are  also  location  dependent. 
Quantitative  assessment  of  these  processes  is,  to  say  the  least,  lacking. 

There  are  still  some  distances  from  a  full  understanding  of  these  processes 
by  the  engineering  community,  particularly  in  the  development  of  analytical 
models  which  describe  and  characterize  them  in  a  way  that  can  be  used  in 
design.  The  overwhelming  difficulties  do  not  appear  to  be  technical  but 
rather  in  overcoming  preconceived  ideas  and  theories  that  do  not  lend  them¬ 
selves  to  a  consistent  treatment  of  the  physical  phenomena  in  high  speed  impact. 

One  of  the  weaknesses  of  classical  physics  and  mechanics  when  used  to 
characterize  the  behavior  of  solid,  liquid  and  gas  is  that  the  theoretical 
formulation  relies  upon  a  knowledge  of  the  constitutive  relation.  As  a  result, 
a  hierarchy  of  theories  has  been  developed,  each  is  being  applicable  over 
a  certain  range  of  stress,  strain,  strain  rate,  temperature  and  material  type 
such  as  elastic  or  plastic  solid,  inviscid  or  viscous  fluid,  perfect  or 
Van  der  Waal  gas,  etc.  If  in  a  single  event,  material  behavior  includes  two 
or  more  of  the  aforementioned  types,  then  the  adjoining  of  the  different 
theories  and/or  solutions  in  itself  becomes  a  monumental  task.  There  is  no 
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hope  for  uniqueness  unless  a  unified  theory  can  be  developed  that  can  account 
for  all  material  behavior  without  making  an  a  priori  assumption  on  the  knowl¬ 
edge  of  the  constitutive  relation. 

Despite  the  extensive  research  efforts  made  on  hypervelocity  impact,  it 
is  not  difficult  to  include  from  past  studies  [4-6]  that  the  traditional 
classical  physics  and  mechanics  approaches  have  not  led  to  any  in-depth  under¬ 
standing  of  the  impact  phenomenon.  Recognized  is  the  fact  that  a  deformable 
solid  can  undergo  several  stages  of  extreme  mechanical  behavior  when  subjected 
to  hypervelocity  impact.  In  the  vicinity  of  initial  contact,  the  solid  behaves 
very  much  like  a  fluid  because  of  the  very  high  elevation  of  local  pressure. 

It  is  tempting  at  this  stage  to  model  the  material  behavior  by  Stokesian  fluid 
or  to  use  the  Navier-Stokes  theory  of  viscous  flow.  In  so  doing,  however,  the 
fluid/solid  interaction  and/or  coupling  effects  are  left  out  as  the  transforma¬ 
tion  of  solid  to  fluid  occurs  only  lodatty  while  the  surrounding  projectile/ 
target  material  is  still  in  the  solid  state.  The  inclusion  of  material  rigidity 
effects  calls  for  the  application  of  other  theories,  say  the  Prandtl -Reuss 
theory  of  elastic-plastic  flow  in  conjunction  with  the  von  Mises  yield  condi- 
tion  .  Additional  assumption  on  fracture  criteria  may  have  to  be  made  as 
failure  of  the  projectile/target  system  in  the  final  stages  can  involve  crack 
propagation.  Since  the  self-contained  theories  for  materials  at  different 

★ 

Such  an  approach  [1]  has  been  developed  in  recent  years  and  applied  to 
explain  a  number  of  previously  not  understood  physical  phenomena  in  material 
damage.  Predictions  have  been  made  for  the  SAFC-40R  steel  [2]  and  6061 -T6 
aluminum  [3]  and  the  results  agreed  very  well  with  experiments. 

** 

Generalization  of  this  condition  to  include  dilatational  effects  would  lead  to 
inconsistency  in  the  plasticity  theory  that  assumes  the  uniaxial  data  to 
coincide  with  the  effective  stress  and  effective  strain  curve.  One  of  the 
inherent  shortcomings  of  all  classical  mechanics  theories  involving  non¬ 
linearity  and/or  dissipation  is  their  inability  to  translate  uniaxial  data  to 
multiaxial  states  without  invoking  serious  physical  limitations. 
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stages,  yield  conditions  and  fracture  criteria  are  unrelated  ,  their  end 
results  would  be  independent.  Any  attempts  made  to  connect  the  individual  solu¬ 
tions  will  involve  arbitrariness  and  conditions  that  may  not  be  consistent  with 
physical  observations.  In  fact,  a  three-phase  equation  of  state  for  metals 
has  been  proposed  [3]  to  account  for  material  in  the  solid-liquid  and  liquid- 
vapor  region.  Use  was  made  of  thermodynamics  and  many  assumptions  that  are 
subject  to  experimental  verification.  The  work  stopped  short  at  predicting  the 
hydrodynamic  behavior  of  a  real  material  and  called  for  the  incorporation  of 
mechanical  constitutive  relations.  To  reiterate,  material  behavior  under  hyper¬ 
velocity  impact  cannot  be  consistently  explained  by  the  classical  piecemeal 
approach . 

Presented  in  this  work  is  a  methodology  that  can  model  the  extremes  of 

material  behavior  in  hypervelocity  impact.  Constitutive  relations  for  each 

element  in  the  projectile/target  system  are  developed  according  to  the  rate  of 

energy  transmission  and  dissipation.  All  forms  of  energy  involving  heat,  sound, 

phase  transformation,  etc.  are  included  in  the  theory  without  neglecting  their 

★  ★ 

coupling  effects.  Thermodynamics  and  heat  transfer  should  be  an  inherent 
part  of  mechanical  deformation.  Only  in  this  way,  a  complete  and  consistent 
description  of  the  material  in  the  solid,  liquid  and  gaseous  state  can  be  given. 
This  knowledge  is  pertinent  for  developing  subscaling  relations  for  the  hyper¬ 
velocity  impact  of  projectile/target  systems.  Compared  are  the  damage  behavior 

★ 

Energy  used  and  dissipated  during  phase  transformation  of  material  in  solid 
state  to  the  fluid  state  for  example  would  not  be  accounted  for. 

★  ★ 

The  reason  that  thermal /mechanical  coupling  effects  have  not  been  satisfac¬ 
torily  treated  up  to  now  is  because  thermodynamics,  heat  transfer  and  mechan¬ 
ical  deformation  are  regarded  as  separate  individual  disciplines  in  classical 
physics  and  mechanics. 
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of  full  scale,  half  scale  and  quarter  scale  projectile/target  systems.  The 
variations  of  geometric  and  kinematic  parameters  with  target  damage  are  also 
discussed  as  the  solid  transforms  into  fluid  in  a  local  region.  As  the 
complete  time  ‘history  of  material  behavior  under  hypervelocity  impact  is  made 
available,  sub-scaling  tests  can  then  be  designed  to  simulate  the  conditions 
of  the  full  scale  projectile/target  system. 

II.  ENERGY  DENSITY  THEORY 

Energy  can  be  converted  from  one  form  to  another  and  this  process  is  in¬ 
herently  irreversible.  The  conversion  of  mechanical  work  to  heat,  for  example 
must  necessarily  involve  energy  dissipation  that  is  unrecoverable.  Even  if 
the  process  were  carried  out  adiabatically  by  insulating  the  system,  local 
irreversibility  can  still  prevail  as  the  collision  of  molecules  will  involve 
energy  loss.  That  is,  each  element  within  the  system  can  be  disturbed  per¬ 
manently.  Local  and/or  global  temperature  change  is  the  result.  The  a  priori 
assumption  that  the  different  forms  of  energy  such  as  heat,  sound,  etc.  can  be 
summed  will  encounter  fundamental  difficulties  when  analyzing  coupling  effects 
Nonlinearity  in  its  entirety  cannot  be  approached  by  adding  more  terms. 

2.1  Energy  Conversion 

Without  loss  in  generality,  the  energy  density  theory  [1,2]  considers 

energy  in  two  basic  forms.  The  portion  which  has  already  dissipated  and  the 

remainder  that  is  still  available.  Energy  conversion  involves  only  the 
★ 

exchange  of  surface  energy  density  (dW/dA).  (i  =  C,n,c)  and  volume  energy 

* 

As  a  special  case,  consider  the  equilibrium  of  a  spherical  liquid  bubble  of 
radius  R  under  uniform  pressure  p.  Equation  (1)  yields  the  well  known  result 
pR  =  2-f  since  dW/dV  =  p,  dV/dA  =  R/2  and  dW/dA  =  y,  the  surface  tension  of 
the  liquid. 


density  dVi/dV  at  a  given  time  t,  i.e., 


,dW,  ,d\/v  dW 

'dA'i  ^  dA J  i  d V  * 


i  =  f. ,  n  , t 


(1) 


Here,  £,  ,n  and  r,  refer  to  the  plane  of  homogeneity  and  are  different  from  the 
reference  axes  x ,y  and  z.  On  this  plane,  the  surface  energy  density  is  assumed 
to  be  the  same  in  all  three  mutually  perpendicular  directions  for  a  homogeneous 
and  isotropic  sol  id: 


/  dWv  _  /  d It ,  _  ,  d W \  <  d W \ 

dAJ f,  '  “  'dA' c  '  Wo 


(2) 


The  value  (dW/dA)Q  is  referred  to  the  uniaxial  case.  The  condition  in  equation 
(2)  as  illustrated  in  Figure  1,  of  course,  must  be  modified  for  materials  whose 
internal  structure  may  possess  preferred  directions.  What  is  important  is  a 
knowledge  of  the  relation  between  (dW/dA).  (i  =  T.n.t)  that  establishes  the 
plane  of  homogeneity.  The  quantity  dW/dV  in  equation  (1)  is  a  scalar  being 
invariant  under  rotation  of  coordinate  axes.  It  follows  that 


(*1)  =  (dVv  =  (dV, 

'dA'1-;  'dA'  n 


dA 'c  'dA  o 


The  change  of  volume  with  surface  (dV/dA).  (i  =  s.o.c)  on  the  plane  of  homo¬ 
geneity  are  also  connected.  Equation  (3)  determines  the  relation  between 
(C»n,c)  and  (x,y,z)  because  the  strain  components  in  (dV/dA ) ^  will  be  known 
from  a  given  deformation  field.  Refer  to  [2]  for  details. 

2.2  Strain  Function 

The  energy  density  theory  requires  only  a  knowledge  of  dV/dA  as  a 
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Figure  1.  Plane  of  Homogeneity:  Equal  Surface  Energy 

Density  for  an  Isotropic  and  Homogeneous  Body 


function  of  uniaxial  strain  for  the  solid,  liquid  or  gas  at  the  initial 
reference  state.  For  the  fluid,  strain  may  be  expressed  in  terms  of  the 
ratio  of  volume  change  to  volume,  i.e.,  AV/V.  The  subsequent  behavior  of 
the  system  is  determined  directly  from  the  load-time  history.  The  strain 
function 


f(0  =  I  x  3 a  dc 


plays  a  fundamental  role  in  the  energy  density  theory  from  which  the  volume- 
energy  function  in  equation  (1)  can  be  obtained  from  a  uniaxial  strain  stat 


The  parameter  x  can  be  adjusted  such  that  f(e)  may  be  identified  with  the  uni¬ 
axial  stress  o(e).  This,  however,  is  not  a  necessary  requirement  in  the  energy 
density  theory.  In  the  case  of  multiaxial  strain,  equations  (4)  and  (5)  may  be 
referred  to  the  plane  of  homogeneity  and  combined  to  yield 


dW  =  rf  wdV 


>!  ideidei  =  e»h  or  t  with  no  sum  on  i) 


Because  of  equation  (3),  it  suffices  to  compute  dW/dV  for  i=^  as  follows: 


dW  _  rr  ,dV  j 
dV  ^  dA^C  £ 


To  be  emphasized  is  that  even  though  equation  (7)  contains  only  the  normal 
strain  dW/dV  is  computed  for  an  element  in  a  multiaxial  strain  state 
involving  the  influence  of  all  the  strain  components  via  (dV/dA)  .  This  pro¬ 
vides  an  unique  means  of  identifying  uniaxial  and  multiaxial  data  without 

★ 

simplifying  the  physics  of  the  problem  . 

The  strain  function  defined  in  equation  (4)  corresponds  to  an  uniaxial 
strain  state.  In  general,  it  may  depend  on  the  strain  components  e^,  etc. 
For  an  isotropic  and  homogeneous  solid  in  two  dimensions,  the  strain  function 
may  be  given  by 


The  classical  theory  of  plasticity  assumes  that  the  uniaxial  data  coincide 
with  the  effective  stress  and  effective  strain  curve.  The  dilatational  effects 
are  thus  neglected  and  cannot  be  accounted  for  by  adding  them  into  the  yield 
condition  as  it  would  lead  to  inconsistency  in  the  formulation. 


i 

\ 

) 

(8)  i 


The  functions 


F  =  F[(dV/dA)5,€e],  G  =  G[(dV/dA)c,e^]  (9) 

can  be  determined  from  the  equations  of  motion. 

2.3  Thermal/Mechanical  Interaction 

Classical  thermodynamics  associates  irreversibility  with  entropy  and  heat 
dissipation.  Since  energy  dissipation  may  take  place  in  so  many  different 
forms,  it  is  not  justified  to  assume,  in  general,  that  dissipation  is  asso¬ 
ciated  exclusively  with  heat  loss.  Referring  to  Figure  2,  the  area  under  opq 
is  dW/dV  or  W  and  is  divided  into 

V  =  ((j^p  *  A  =  ^  1 0 

such  that  V  corresponds  to  opp'  and  A  to  pqp'.  The  corresponding  strains  are 

e  =  tp  +  t  (11) 

the  dissipation  energy  density  V  must  be  a  positive  definite1  function. 


+This  condition  can  be  used  as  a  guideline  for  selecting  the  appropriate  grid 
size  and  time  increment  in  numerical  analysis. 
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Figure  2.  Uniaxial  Stress  and  Strain  Response 


Because  of  conservation  of  energy,  it  follows  that 


dW  =  dD  +  dA, 


dO  n 
dt 


02) 


The  rate  of  change  of  V  with  time  is  required  to  be  a  monotonical ly  increasing 
function  of  time. 

★ 

Now,  let  A  depend  on  the  function  H  and  uniaxial  strain  e  such  that 


A  =  A  (H,0 


(13) 


For  a  multiaxial  strain  state,  c  on  the  plane  of  homogeneity  may  be  used 
instead  of  e.  ^ 
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Applying  the  rules  of  differentiation,  equation  (13)  gives 


dA  .  (§£)*,  ♦  (|A)Hd£ 


Defining 


0  =  ■> 


- 

0  "  3  e  1/ 


and  differentiating  equation  (15),  it  is  found  that 


3c  H  Ud H'e 


With  the  aid  of  equations  (12),  (14)  and  (15),  a  definition  of  the  H  function 


is  obtained: 


-dH  .  f 


The  negative  sign  arises  as  work  is  assumed  to  be  done  on  the  system  rather 


than  as  done  by  the  system  in  classical  thermodynamics.  Equation  (17)  resem¬ 


bles  the  classical  definition  of  entropy  as  related  to  heat  Q  and  temperature 


T.  The  quantity  o  may  be  scaled  in  °K  as  in  the  case  of  T  and  it  becomes  T 


only  when  energy  dissipation  occurs  all  in  the  form  of  heat,  i.e.,  V  =  Q. 


For  a  finite  increment  of  strain  Ac,  it  can  be  shown  from  equations  (16)  and 


(17)  that 


For  an  element  in  a  multiaxial  stress  or  strain  state,  equation  (18)  should  be 
referred  to  the  plane  of  homogeneity  so  that 


A0 

0 


aV 


(19) 


The  temperature  change  at  every  point  in  the  system  can  thus  be  found  without 

★ 

invoking  the  separate  discipline  of  heat  transfer  .  Equation  (19)  reveals 
the  time  nature  of  thermal /mechanical  coupling  that  inherently  involves  energy 
dissipation. 


2.4  Phase  Change 

If  the  rate  of  energy  transfer  in  a  unit  volume  of  material  becomes 
sufficiently  intense,  phase  transitions  in  the  form  of  melting,  vaporization 
and  sublimation  can  take  place.  The  transition  is  strain  rate  dependent,  an 
effect  that  is  not  considered  explicitly  in  classical  thermodynamics.  Since 
transition  would  alter  the  order  and  disorder  of  the  system,  it  would  appear 
as  an  inflection  point  on  the  H  versus  time  curve.  A  useful  quantity  for  the 
quantitative  assessment  of  phase  change  is  the  strain  rate  of  dissipation 
energy  density  defined  as 


aV  aQ 
Ag  "■  AV 


(20) 


which  reduces  to  aQ/aV  only  as  a  special  case.  The  values  of  aQ/aV  correspond¬ 
ing  to  the  onset  of  phase  change  for  metals  can  be  found  in  handbooks. 


Heat  transfer  by  radiation,  convection  and  conduction  will  all  be  included  in 
V.  There  is  no  need  to  invoke  assumptions  on  the  way  with  which  heat  and 
temperature  or  temperature  gradient  are  connected. 
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III.  PROBLEM  STATEMENT  AND  RESULTS 


Suppose  that  a  tungsten  projectile  traveling  at  9,000  m/sec  or  30,000 
ft/sec  comes  in  contact  with  a  5083  aluminum  target  plate  as  shown  in  Figure 
3.  Depending  on  the  time  interval  at  which  the  physical  event  is  being 
reviewed,  the  material  will  undergo  deformation,  damage,  phase  transition, 
etc.  The  application  of  the  energy  density  theory  involves  finding  the  dis- 


Center 


Figure  3.  Hyperveloci ty  Impact  of  Projectile/Target  System 

placements  that  yield  the  strain  field  from  which  the  plane  of  homogeneity 
can  be  located.  Unlike  classical  continuum  mechanics,  the  tractions  and 
inertia  forces  will  interact  as  the  rate  change  of  volume  with  surface  is  no 
longer  assumed  to  vanish,  i.e., 


(21) 


?  ,  32ui  ,dV, 

Ti  =  °ijnj  +  p  IF"  Wn 


in  which  p  is  the  mass  density  and  n  the  normal  to  the  plane  on  which  the 
n 

tractions  T.  prevail.  Here,  o..,are  the  stress  components  that  can  be  related  to 

•  I J 

the  strain  function  for  each  strain  state  through  a.  The  displacement  components 
are  denoted  by  u^ .  It  is  important  to  recognize  that  even  though  (dV/dA)^  in 
equation  (21)  may  be  small,  the  term  p  a2u^/at2  may  be  large  in  dynamic  problems 

3.1  Initial  Material  Response 

As  remarked  earlier,  only  a  knowledge  of  the  initial  mechanical  and 

thermal  properties  is  required.  They  are  outlined  in  Table  1  for  the  tungsten 

projectile  and  5083  aluminum  target.  The  quantities  in  Table  1  will  be 

defined  as  the  discussion  proceeds. 

TABLE  1.  Reference  Mechanical  and  Thermal  Properties 
of  Projectile  and  Target 


Parameter 

Tungsten 

5083  Aluminum 

E(GPa) 

345 

64.2 

oys(MPa) 

1 ,062 

132.5 

o^(MPa) 

1  ,347 

260.8 

(dW/dV)c(MPa) 

15.03 

9.89 

( dW/dA ) c( N/m) 

6.18xl04 

2.09xl05 

p(kg/m3) 

19,350 

2,700 

C  ( J/kg°K) 

145.7 

900 

T  (°K) 

3,683 

1  ,200 

m 

aQ( J/kg) 

4.93x105 

8 .10x1 05 

AQ/AV(Pa) 

9.64xl07 

2 . 2 1  x 1 07 

For  the  sake  of  illustration,  Figures  4  and  5  display  the  family  of 
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Stress  o  x  10  (MPa) 


stress  and  strain  curves  for  tungsten  and  5083  aluminum  as  the  strain  rates 
are  varied.  They  are  referred  to  as  the  material  data  bank  such  that  the 
base  material  serves  as  the  reference  state.  Refinement  of  the  data  can  be 
achieved  by  introducing  a  floating  bank  to  the  desired  degree  of  accuracy. 

The  idea  is  that  each  point  (a  ,e.)  on  the  plane  of  homogeneity  can  be 
uniquely  and  accurately  identified  with  a  uniaxial  stress  and  strain  state 
or  [ ( dV/dA)  ,c  ] .  The  actual  response  of  a  given  element  in  the  projectile/ 
target  system  may  not  coincide  with  any  of  the  curves  shown  in  Figures  4  and  5 
It  will  most  likely  follow  a  path  that  intersects  the  curves  at  several  points 
because  the  probability  of  any  elements  deforming  at  a  constant  strain  rate 
under  impact  is  small.  As  it  will  be  shown  subsequently,  the  stress  and 
strain  curves  for  elements  near  the  region  of  impact  will  in  no  way  resemble 
those  in  Figures  4  and  5  for  tungsten  and  aluminum  in  the  solid  state.  In 
fact,  they  will  behave  as  viscous  fluid  because  local  melting  has  occurred. 
This  is  illustrated  in  Figure  6  where  the  initial  responses  of  the  projectile 
and  target  material  are  covered  by  the  family  of  linear  curves.  The  dotted 
line  is  the  actual  response  of  an  element  in  the  fluid  state  determined  auto¬ 
matically.  Phase  transition  is  accounted  for  by  the  energy  density  theory. 
There  is  no  need  to  guess  or  anticipate  the  form  of  the  constitutive  relations 
as  required  in  the  classical  approach.  The  shapes  of  the  stress  and  strain 
curves  selected  in  the  material  data  bank  are  not  essential  as  long  as  all  the 
points  in  the  versus  domain  are  uniquely  identified  with  [ (dV/dA) ^ ,  e.]. 
The  objective  is  to  mathematically  establish  an  one-to-one  correspondence 
between  the  stress  and/or  strain  state  of  each  element  and  its  equivalent 

★ 

Each  curve  refers  to  a  different  strain  rate. 
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Figure  6.  Response  of  Viscous  Fluid  Determined 
From  Material  Data  Bank  For  Solid 


uniaxial  state.  It  is  simply  a  bookkeeping  procedure.  There  is  no  longer  the 

★ 

need  to  obtain  uniaxial  data  for  different  strain  rates  except  for  the 
initial  response  to  be  used  as  a  reference.  The  enetgy  demlty  theoty  can 
analytically  detexinine  the  complete  A-tte-ii  and  Attain  behavior  iox  any 
Attain  tatCA  ^tom  only  a  knowledge  ofi  the  initial  Young '6  moduli- 4  and  the 

■specified  loading  note.  This  has  been  done  for  the  6061  -T6  aluminum' 

-3  -4  -1 

stretched  uniaxially  at  a  strain  rate  of  10  and  10  sec  and  the  predic¬ 
tion  is  within  5%  of  the  experimental  data  [8]. 


3.2  Image  Mapping 

The  projectile/target  system  is  discretized  by  application  of  the  twelve 
(12)  node  isoparametric  finite  elements.  Since  the  geometry  in  Figure  1 


For  strain  rates  of  10  sec"  and  higher  considerable  difficulties  are 
encountered  experimentally  in  uniaxial  tests.  A  review  of  the  subject  can 
be  found  in  [7]. 
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possesses  one-half  symmetry,  it  suffices  to  use  20  elements  and  138  nodes 
as  indicated  in  Figure  7(a).  As  the  boundary  conditions  along  the  plane  of 
symmetry  cannot  always  be  adequately  satisfied  numerically,  the  method  of  image 
mapping  will  be  adopted.  Sixteen  additional  elements  are  added  that  simply 
reflects  the  displacement  field  symmetrically  across  the  center  plane, 

Figure  7(a).  This  will  obviously  provide  more  accurate  results  for  elements 
no.  15  and  17.  Particular  attention  will  be  given  to  the  results  at  points 
1,  2  and  3  in  the  target  and  4,  5  and  6  in  the  projectile,  Figure  7(b).  All 
quantities  referred  to  these  locations  will  be  averaged  over  a  circular  region 
with  diameter  0.02  cm.  Computed  are  dW/dV  and  dW/dA  contours  in  the  pro¬ 
jectile/target  system  for  different  time  increments.  Figure  8  displays  the 
volume  energy  density  contours  at  t  =  2  nsec.  The  contours  are  seen  to  be 
densely  packed  near  the  corner  with  increasing  magnitude.  This  is  not 
surprising  because  reentrant  corners  tend  to  concentrate  energy.  Even  at  this 
very  early  stage  of  contact,  energy  has  transferred  to  all  elements  in  the 
target.  There  prevails  a  five  order  of  magnitude  difference  in  the  intensity 
of  the  twenty  one  (21)  dW/dV  contours  displayed  in  Figure  8.  Similarly,  the 
surface  energy  density  contours  can  also  be  plotted  as  given  in  Figure  9  at 
t  =  2  nsec.  Their  distribution  is  similar  to  that  of  dW/dV  in  Figure  8.  The 
spread  in  the  difference  of  the  intensity  of  dW/dA,  however,  has  increased  by 
one  order  of  magnitude,  i.e.,  the  value  of  dW/dA  for  contour  no.  1  differed 
from  that  of  no.  21  by  six  order  of  magnitude.  The  surface  energy  density  is 
expected  to  dominate  at  initial  impact  because  energy  is  transmitted  across 
the  contact  surface  between  the  project  and  target  before  it  can  spread  vol¬ 
ume  wise.  An  enlarged  view  of  the  dW/dA  contours  near  the  corner  is  shown 
in  Figure  10.  This  provides  an  insight  into  the  shape  of  local  region  within 


(a)  Nurobe 


Figure  9.  Constant  Surface  Energy  Density 
dW/dA  Contours  at  t  =  2  nsec 
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which  the  material  may  change  phase.  As  time  elapses  to  t  =  4  nsec,  energy 
is  being  continuously  transmitted  from  the  projectile  into  the  target.  This 
is  evidenced  by  the  increase  in  the  intensity  of  the  dW/dV  and  dW/dA  contours 
exhibited  in  Figures  11  and  12,  respectively.  An  one  order  of  magnitude 
increase  in  the  volume  and  surface  energy  density  can  be  observed.  Concen¬ 
tration  of  energy  remains  in  the  vicinity  of  the  corner.  Figure  13  gives  the 
details  of  the  dW/dA  contours  for  t  =  4  nsec. 

IV.  SUB-SCALING  OF  PROJECTILE/TARGET  SYSTEM 

Prior  to  performing  sub-scaling  tests  of  projectile/target  systems,  it  is 
necessary  to  establish  the  conditions  or  criteria  under  which  geometric, 
kinematic,  material,  damage  and  time  parameters  are  to  be  proportioned  or 
scaled.  Because  energy  dissipation  during  hypervelocity  impact  plays  a 
dominant  role,  scaling  of  the  aforementioned  parameters  will  be  highly  non¬ 
linear.  For  the  purpose  of  gaining  insight  into  the  physical  problem,  three 
different  projectile/target  systems  will  be  analyzed. 

4.1  Geometric  Proportionality 

Shown  in  Figures  14(a),  14(b)  and  14(c)  are  three  different  projectile/ 
target  systems  which  shall  be  referred  to  as  full  scale  (1:1),  half  scale  (1:2) 
and  quarter  scale  (1:4),  respectively.  Geometric  proportionality  is  assumed 
such  that 


A_a_a  C_c  y 

B  "  b  6  ’  0  '  d  "  6 


(22) 


The  masses  m  ,  in,  , 
p  h 


triq  and  velocities 


V 


V 


q 


tan  alter  depending  on  the 
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Contour  Values  of  dW/dV(Pa) 
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Figure  11.  Constant  Volume  Energy  Density 
dW/dV  At  t  =  4  nsec 


Contour  Values  of  dW/dA(Pa-ro) 
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Figure  12.  Constant  Surface  Energy  Density 
dW/dA  Contours  At  t  =  4  nsec 
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Contour  Values  of  dW/dA(Pa-m) 
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Figure  13.  Constant  Surface  Energy  Density  dVI/dA 
Contours  On  Enlarged  Corner  Region  Of 
Projectile/Target  System  At  t  =  4  nsec 
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(a)  Full  scale  (1:1)  (b)  Half  scale  (1:2)  (c)  Quarter  Scale  (1:4) 

Figure  14.  Sub-scaling  of  Projectile  Target  System 


extent  with  which  kinetic  energy  and/or  momentum  change  affect  material 
damage . 

Suppose  that  the  systems  in  Figures  14(a)  to  14(c)  are  scaled  to  preserve 
the  kinetic  energy  in  addition  to  satisfying  the  conditions  in  equation  (22). 
This  then  requires 


Ivi 


\  m  v* 
2  q  q 


(23) 


for  the  initial  time  step.  Solving  for  the  velocities  of  the  sub-scale  system 
in  terms  of  the  full  scale  projectile  velocity,  it  is  found  that 
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m  0.5 

V,  =  (-^)  v 

h  mh  P 


m  0.5 

v  =  (— ^)  v 
Q  iy  P 


On  the  other  hand,  if  momentum  is  to  be  conserved  then 


(24) 


m  v  =  m.  v„ 
P  P  bp 


m  v 

q  q 


(25) 


This  yields 


m  1 

v.  =  (-£•)  v 
h  p 


m  1 

v  -  l  -£-)  V 
q  'ni  ;  f 


(26) 


In  practice,  the  exponent  on  the  mass  ratio  will  neither  be  0.5  or  1.0.  That  is 
geometric  proportionality  will  not  conserve  the  projectile  kinetic  energy  or 
momentum.  At  any  given  instance  of  time,  there  will  be  a  two  parameter 
scaling  relation  of  the  form 


m_  nl 


v' 


(27) 


for  the  half-scale  system  and 


m 

■ 

q 


(28) 


for  the  quarter-scale  system.  The  quantities  k-j  ,  k^  and  will  depend 

on  the  damage  in  the  projectile/target  system  that  changes  with  time. 

4  .2  Response  in  Target 

If  the  three  systems  in  figures  14(a)  to  14(c)  are  to  experience  approxi- 
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mately  the  same  damage  at  a  given  time,  then  the  velocities  vh  and  vq  must 
be  scaled  accordingly.  After  several  trial  runs,  the  initial  velocities  for 
the  half  scale  and  quarter  scale  model  in  Table  2  were  chosen.  They  are 
lower  than  that  of  the  full  scale  system.  Otherwise,  the  scaled  down 
target  will  experience  more  damage  at  a  given  time.  For  the  -same  ma.te'ccaO, 
the  energy  concentration  will  be  more  severe  for  the  smaller  projectile/ 
target  system  because  the  distance  between  the  corners  are  reduced.  Figures 
15,  16  and  17  show  the  variations  of  compressive  stress  with  compressive 

TABLE  2.  Projectile  Velocity  at  t  =  4  nsec 
for  Different  Scaled  Models 


Initial 

Velocity  at 

Scale 

Velocity 

t  =  4  nsec 

(m/sec) 

(m/sec) 

Full  (1:1) 

9,000 

8,526.7 

Half  (1:2) 

4,000 

3,836.5 

Quarter  (1:4) 

1  ,800 

1  ,690.0 

strain  c  for  the  full  scale,  half  scale  and  quarter  scale  system  at  t  =  4 
nsec  at  points  1,  2  and  3  as  indicated  in  Figure  7(b),  respectively.  The 
response  corresponds  to  that  for  viscous  fluid.  Point  2  being  directly  under 
the  periphery  of  the  projectile  experiences  higher  stress  and  strain.  Next 
in  line  is  point  3  which  is  near  the  free  surface  of  the  target.  The  term¬ 
inal  points  of  the  a  versus  ef  curves  in  Figures  16  and  17  can  be  elevated 
closer  to  those  in  Figure  15  for  the  full  scale  system  if  the  initial  veloc¬ 
ities  v,  and  v  in  Table  2  are  increased  accordingly.  Only  the  initial 
h  q 

response  of  the  base  material  were  used  to  obtain  the  results  in  Tigures  15 
to  17.  The  corresponding  numerical  values  are  given  in  Table  3  for  t  =  2 
and  4  nsec  . 


Pefer  to  detailed  discussion  in  Section  V. 
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Compressive  stress  a  xlO  (Pa) 


Compressive  stress 
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Figure  16.  Stress  and  Strain  Response  At  Point  1,  2 
And  3  In  Target  For  t  =  4  nsec  And 
Half  scale  (1:2) 
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Figure  17.  Stress  And  Strain  Response  At  Point  1,  2 
And  3  In  Target  For  t  =  4  nsec  And 
Quarter  Scale  (1 :4) 
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TABLE  3.  Stress  and  Strain  on  Plane  of  Homogeneity 
at  Point  1 ,  2  and  3  in  Target 


Point  1  _ Point  2  Point  3 


_3 

Time  e^xlO 

t(nsec)  (cm/cm) 

o^xlO8 

(Pa) 

e^xlO'3 
(cm/ cm) 

o^xlO8 

(Pa) 

^xlO"3 

(cm/cm) 

1 

o^xl  0 
.  (Pa) 

Full  Scale  (1  :1 ) 

2  -0.08965 

-0.08745 

-0.1758 

-0.1829 

-0.1134 

-0.1303 

4  -0.1824 

-0.3056 

-0.3533 

-0.5994 

-0.2117 

-0.4324 

Half  Scale  (1:2) 

2  -0.07969 

-0.07813 

-0.1563 

-0.1634 

-0.1008 

-0.1164 

4  -0.1622 

-0.2745 

-0.3142 

-0.5384 

-0.1882 

-0.3883 

Quarter  Scale  (1:4) 

2  -0.07172 

-0.07095 

-0.1407 

-0.09069 

-0.1057 

4  -0.1460 

-0.2827 

-0.3559 

By  determining  the  path  of  local  unloading,  the  available  energy  density 

can  be  distinguished  from  the  dissipation  energy  density  for  each  load  step. 

The  quantity  V  as  defined  in  equation  (10)  can  thus  be  found.  Keeping  in  mind 

★ 

that  V  and  dP/dt  must  be  positive  .  The  correct  feature  and  trend  of  V  is  shown 
in  Figures  18(a)  to  18(c)  for  the  three  different  systems.  Point  2  dissipates 
more  energy  at  t  =  4  nsec  than  points  1  and  3.  This  trend  remains  valid  for 
t  =  2  nsec  as  given  in  Table  4.  Once  the  time  history  of  P  is  known,  &P  in 
equation  (19)  follows  immediately.  With  the  aid  of  the  data  in  Table  3  for 
and  e^,  the  temperature  o  in  equation  (18)  can  be  calculated.  The  results 
are  plotted  as  a  function  of  time  in  Figures  19(a)  to  19(c)  with  o  being 
★ 

These  conditions  may  not  be  satisfied  on  the  first  trial.  Readjustment  of 
the  finite  element  grid  size  compatible  with  the  time  increment  is  freauently 
necessary.  The  optimum  choice  corresponds  to  diminishing  oscillation  of  V 
wi th  time  for  smal 1  t . 
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equal  to  300°K.  Even  though  all  three  points  1,  2  and  3  are  beyond  the 

melting  tenperature  0m  =  1 ,200°K  for  the  5083  aluminum,  the  strain  rate  of 

energy  dissipation  density  AP/Ae  should  be  checked  to  see  whether  the  mater i a  1 

elements  at  points  1,  2  and  3  have  indeed  melted  or  not.  Refer  to  Table  5 

for  the  numerical  values  of  o-o  .  It  is  conceivable  that  AP/Ae  corresponding 

★ 

to  a  unit  volume  of  material  may  pass  the  threshold  of  melting  only 
momentarily  and  then  drops  below  its  critical  value.  According  to  the 
results  in  Table  4,  this  is  not  the  case.  But  rather,  AP/Ae  at  points  2 
and  3  increased  monotonically  with  time  as  it  passes  the  threshold 
aQ/aV  =  22.1x1 O^Pa  as  given  in  Table  1  for  the  5083  aluminum.  Point  1  remains 
below  the  melting  condition  even  though  its  temperature  is  near  the  melting 
point.  The  transition  of  solid  to  fluid  for  a  unit  volume  of  material  is 
guaranteed  when  both  o  and  aP/Ac  are  maintained  beyond  the  critical  level 
within  a  certain  time  period. 

To  be  noted  is  that  at  the  very  early  stage  of  impact,  the  phenomenon  of 
cooling  and  heating  also  occurs  although  its  influence  is  negligible  in  com¬ 
parison  with  the  time  history  of  the  total  impact  process.  This  is  similar 
to  that  observed  in  [2,3]  for  metal  specimens  loaded  uniaxially  at  low  strain 
rates . 

4.3  Response  in  Projectile 

Upon  impact,  the  elements  in  the  projectile  also  experiences  compressive 
stress  and  compressive  strain.  Referring  to  the  plane  of  homogeneity,  the 
o,  versus  tr  relations  for  points  5  and  6  at  t  =  4  nsec  are  plotted  in 
Figures  21,  22  and  23.  They  represent,  respectively,  the  full  scale,  half 
scale  and  quarter  scale  as  shown  in  Figures  14(a),  14(b)  and  14(c).  Again, 

★ 

A  unit  volume  of  material  aV  used  to  calculate  t. P/a.,  corresponds  to  a  local 
region  0.02  cm  in  diameter. 
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TABLE  4.  Dissipation  Energy  Density  and  Strain  Rate 
Dissipation  Energy  Density  at  Point  1,  2 
and  3  in  Target 


Full  Scale  (1  :1 ) 


2 

0.1625 

4 

1  .3143 

Half  Scale  (1:2) 

2 

0.1303 

4 

1  .0575 

Quarter 

Scale  (1 :4) 

Point  1 


AD/AeXlO 


4.12 

19.96 

3.66 

18.05 


0.1083 
.88 


Point  2 


tiO/  Aexl  0 


0.6496 

4.6087 

0.5204 

3.7091 

0.4318 


6.09 

38.20 

5.32 

34.79 


Point  3 


A V  AexlO 


0.3259 
1  .9972 

0.2605 
1  .6048 

0.2154 


TABLE  5.  Temperature  Change  at  point  1,  2  and 
3  in  Target 


Base  material 


Initial  response 


_ i _ i _ i _ i _ 

•0 . 2  -0.4  -0.6  -0.8 

Compressive  strain  e^xlO  (cm/cm) 

Figure  21.  Local  Stress  And  Strain  Response  At  Point  5 
and  6  in  projectile  for  t  =  4nsec  And 
Full  scale  (1:1) 
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Figure  22.  Local  Stress  And  Strain  Response  At  Point  5 
And  6  In  projectile  For  t  =  4  nsec  And 
Half  Scale  (1:2) 
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Base  material 


/  Initial  response 


= - 1 - — i— - - - 1 _ i _ 

-0.2  -0.4  -0.6  -0.8 

-3 

Compressive  strain  c^xlO  (cm/cm) 

23.  Local  Stress  And  Strain  Response  At  Point  5 
And  6  In  Projectile  For  t  =  4  nsec  And 
Quarter  Scale  (1 :4) 


( 


only  the  initial  response  of  tungsten  or  base  material  is  used.  The  corner 
element  in  the  projectile  corresponding  to  point  6  in  Figure  7(b)  is 
stressed  far  beyond  point  5  just  inside  the  projectile.  Both  curves  labelled 
5  and  6  in  Figures  21  to  23  are  slightly  nonlinear.  This  effect  will  become 
more  obvious  when  damage  begins  to  take  place  at  the  later  stages.  For 
details,  reference  can  be  made  to  Table  6  for  the  numerical  values  of  o  and 
c  including  those  at  t  =  2  nsec.  The  dissipation  energy  density  V  as  a 
function  of  time  is  given  in  Figures  24(a)  to  24(c).  In  all  cases,  V  increased 
slowly  at  first  and  then  rises  very  quickly.  This  trend  is  particularly  pro¬ 
nounced  for  point  6  even  though  the  elapsed  time  is  only  4  n  sec.  The 

TABLE  6.  Stress  and  Strain  on  Plane  of  Homogeneity 
at  Point  5  and  6  in  Projectile 


Point 

6 

Point 

5 

T  ime 

c  xlO'3 

a^xlO^ 

v10"3 

(cm/cm) 

o^xl 08 

t(nsec)  (cm/cm) 

(Pa) 

(Pa) 

Full 

Scale  (1 : ) 

2 

-0.0762 

-0.3981 

-0.0198 

-0.1141 

4 

-0.1635 

-1.3548 

-0.0387 

-0.3653 

Half 

Scale  (1  : 2) 

2 

-0.0677 

-0.3557 

-0.0176 

-0.1019 

4 

-0.1454 

-1  .2168 

-0.0344 

-0.3280 

Quarter  Seal e  (1:4) 

2 

-0.06097 

-0.3230 

-0.01585 

-0.09256 

4 

-0.1309 

-1 .1159 

-0.03097 

-0.3007 

dissipation  at 

point  5  is  relatively  low  as 

shown  by  the 

values  in  Table  7  . 

Plotted  in  Figures  25(a)  and  25(c) 

are  the 

time  histories 

of  the  temperature 

i 

i 

f 

4 


J 


I 


I 
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Figure  24.  Time  History  Of  Dissipated  Energy  Density  At  Point 
5  And  6  In  Tungsten  Projectile 


Figure  26.  Time  History  Of  Strain  Rate  Dissipati 
At  Point  5  And  6  In  Tungsten  Projecti 


TABLE  7.  Dissipation  Energy  Density  and  Strain  Rate 
Dissipation  Energy  Density  at  Point  5  and 
6  in  Projectile 

_ Point  6 _  _ Point  5 _ 

Time  PxlO3  AP/AexlO6  PxlO3  AP/AexlO6 

t(nsec) _ (Pa) _ (Pa) _ ( Pa) _ ( Pa) 

Full  Scale  (1:1) 


2 

0.5151 

12.82 

0.04532 

3.27 

4 

4.0931 

76.70 

0.2756 

24.43 

Half  Scale 

(1:21 

2 

0.41 31 

11.56 

0.03625 

3.30 

4 

3.2967 

69.49 

0.2216 

22.32 

Quarter  Scale  (1 :4) 

2 

0.3435 

10.91 

0.03003 

3.10 

4 

2.7646 

62.30 

0.1855 

20.68 

change  at  points  5  and  6  in  the  projectile.  A  sign  change  in  0-0o  is  again 
observed  for  small  time.  The  temperature  rose  quickly  in  all  the  cases  and 
then  began  to  level  off  for  t  >  2  nsec.  The  plateaus  of  all  the  curves  is 
substantially  below  the  melting  point  of  tungsten  which  occurs  at  3,683°K  as 
given  in  Table  1.  Refer  to  Table  8  for  numerical  values  of  0-0  .  Finally, 
Figures  26(a)  to  26(c)  show  the  variations  of  aP/Ae  with  time.  All  the  curves 
are  below  the  threshold  level  of  aQ/aV  =  96.4xl06Pa  for  tungsten.  No  melting 
takes  place  in  the  projectile.  Table  7  reveals  that  the  highest  aP/ae  value 
is  76.70xl0^Pa  and  this  occurred  at  point  6  for  t  =  4  nsec. 


TABLE  8.  Temperature  Change  at  Point  5  and  6  in 
Projectile 


Point  6  Point  5 

Time  o-o  G-0 

o  o 


2  1,149  979 

4  1  ,059  1  ,007 


4.4  Local  Damage  in  Target 

From  the  results  in  Sections  4.2  and  4.3  it  can  be  concluded  that  only 
the  target  material  next  to  the  projectile  corner  is  being  damaged.  The 
material  in  fact  has  locally  undergone  a  phase  transition  where  the  solid 
has  melted  and  turned  into  fluid.  Such  a  behavior  is  reflected  by  the  stress 
and  strain  curves  in  Figures  15  to  17.  The  shaded  regions  in  Figures  27(a) 
to  27(c)  are  zones  in  which  the  solid  has  melted.  Their  volumes  and  areas 
denoted  by  Vm  and  Am>  respectively,  are  estimated  and  given  in  Table  9.  The 
ratio  Vm/ Am  is  the  largest  for  the  full  scale  system  and  the  smallest  for  the 
quarter  scale  system.  Their  relative  difference  will  change  with  time  in  a 
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Molten  Local  Region  In  Tarqet/Projectile  System  For 
T i me  Different  Scale  Size  At  t  =  4  nsec 


TABLE 

9.  Size 
t  =  4 

of  Melted  Target  Region  at 
nsec 

Volume 

Area 

Ratio 

Impact  Veloci t 

Scale 

Vc"|3> 

Am(cm2) 

V  /A  (cm) 
m  in  ' 

m/ sec 

Full  (1:1) 

0.1325 

0.9843 

0.1346 

9,000 

Half  (1:2) 

0.0231 

0.4528 

0.0510 

4,000 

Quarter  (1:4) 

0.00391 

0.1870 

0.0209 

1 ,000 

nonlinear  fashion.  The  complexities  arising  from  the  interaction  of  the 
geometric  and  kinematic  parameters  with  damage  have  been  clearly  established. 
Any  conclucions  prior  to  a  complete  knowledge  of  the  time  damage  process  of 
the  different  scaled  down  systems  would  be  premature. 

V.  ADDITIONAL  REMARKS  ON  SCALING 

The  foregoing  results  provide  sufficient  evidence  to  the  fact  that  damay 
caused  by  hypervelocity  impact  depends  upon  the  complex  interaction  of 
several  variables.  There  exists  no  simple  laws  of  scaling.  Geometric  propoi 
tlonallty  alone  will  not  preserve  the  damage  behavior  of  the  projectile/ 
target  system.  This  can  be  demonstrated  by  solving  the  half  scale  problem 
with  the  projectile  traveling  at  the  same  velocity  as  in  the  full  scale  case 
Refer  to  the  situations  depicted  in  Figures  28(a)  and  28(b) .  For  the  same 
projectile  velocity  of  9,000  m/sec  or  30,000  ft/sec,  it  will  be  shown  that 
the  target  damage  in  the  half  scale  system  will  in  no  way  resemble  that  of 
the  full  scale  system  at  t  =  4  nsec.  This  difference  will  increase  with  tim 


Referring  to  the  data  in  Table  10  and  the  sketches  in  Figures  28(a)  and 
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Figure  28.  Cylindrical  Projectile  Impacting  Circular  Target  Plate 
At  Hypervelocity  With  Damage  At  t  =  4  nsec 
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TABLE  10.  Damage  Zones  for  Full  and  Half  Scale 
Projectile  With  the  Same  Initial 
Velocity  of  9,000  m/sec  and  t  =  4  nsec 


Scale 

Volume 

V  (cm3) 
m 

Area 

A  (cm^) 
nr 

Ratio 

V  /A  (cm) 
m  mv  ' 

Full  (1:1) 

0.1325 

0.9843 

0.1346 

Half  (1:2) 

0.1438 

0.9940 

0.1461 

28(b),  the  damaged  zone  size  for  the  scaled  down  model  is  larger  than  that  of 
the  full  size  model.  The  data  in  Table  11  also  show  that  the  elements  at 
points  1,  2  and  3  in  the  half  scale  mode  are  all  in  the  fluidic  state  while 


TABLE  11.  Dissipation  Energy  Density  and  Strain  Rate 
Dissipation  Energy  Density  at  Point  1,  2 
and  3  in  Target  for  Half  Scale  System  (1:2) 
with  vh  =  9,000  m/sec 


Point  1 

Point  2 

Point  3 

Time 

t(nsec) 

PxlO3 

(Pa) 

AP/ Aexl 06 
(Pa) 

PxlO3 

(Pa) 

AP/AexlO6 

(Pa) 

PxlO3 

(Pa) 

AP/AcxlO6 

(Pa) 

2 

0.6594 

8.36 

2.6342 

12.20 

1 .3191 

10.89 

4 

5.3501 

40.67 

18.7635 

78.27 

8.1179 

68.96 

point  1  remains  unmelted  in  the  full  scale  model  where  AP/Ae  is  below  the 
critical  value  of  aQ/aV  for  tungsten.  Moreover,  the  strain  rate  of  energy 
dissipation  in  the  half  scale  model  is  more  intense.  At  point  3,  a  simple 
calculation  shows  that  ( aP/ac)^ ^/(aP/ae)^ =  2.1  which  is  a  factor  two 
difference. 

It  is  instructive  to  compare  the  normalized  volume  and  surface  area 
ratio  of  the  damage  zones  for  the  full  scale  and  half  scale  system 
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impacted  at  the  same  velocity*  Making  use  of  the  dimensions  specified  and 
the  data  in  Table  10,  the  following  results  are  obtained: 


where  Vq  and  Aq  are  the  respective  volume  and  area  of  the  target.  The 
normalized  damage  zone  for  the  full  scale  and  half  scale  differed  by  a  factor 
of  approximately  two  at  t  =  4  nsec.  The  difference  will  obviously  alter  with 
time . 


The  velocity  of  the  projectile  alone  does  not  control  the  extent  of 
target  damage.  In  other  words,  the  same  degree  of  damage  can  be  produced  in 
a  sub-scale  model  even  though  the  projectile  may  be  traveling  at  a  much  lower 
speed.  Returning  to  the  example  discussed  earlier  in  Section  IV  where  the 
velocity  of  the  half  scale  projectile  was  4,000  m/sec  rather  than  9,000  m/sec. 
The  damage  zone  at  the  same  instance  is  more  localized  as  shown  in  Figure 
27(b).  Making  use  of  the  data  in  Table  9,  the  normalized  volume  and  area 
ratio  of  the  full  scale  and  half  scale  become 


<Wl:2 

<Wl:l 


=  0.348 


<VA„>1:2 

‘VVlil 


=  0.920 


Although  the  volume  ratio  is  much  less  than  unity,  the  surface  one  ratio  is 
nearly  one.  The  relative  size  of  the  damage  area  is  being  preserved  at 
t  =  4  nsec.  Comparison  of  (aP/ae)^  ^  anc*  ( At>/ Ac ) -i  . ^  can  also  be  made  for 
the  element  say  again  at  point  3.  From  the  data  in  Table  4,  (aP/Ac)^/ 

( ^.D/ zit. ) i . i  =  0.917  at  t  =  4  nsec.  The  strain  rate  of  dissipation  is  being 
preserved.  Nearly  the  same  ratio  also  prevails  for  points  1  and  2.  Should 


the  damage  be  preserved  for  the  entire  time  history,  then  the  time  variable 
must  be  normalized  and  scaled  accordingly. 

No  useful  information  would  be  gained  in  scaling  if  the  extremes  of  the 
geometric  and  kinematic  parameters  were  chosen  while  emphasizing  only  the 
post  mortem  damage  pattern  or  damage  at  any  other  specific  time  instance. 

Beyond  certain  combination  of  the  projectile  speed  and  target  plate  thickness, 
the  resulting  hole  size  in  the  target  will  be  geometrically  proportional  to 
the  projectile  diameter.  This  represents  the  situation  of  overkill  where  the 
target  no  longer  plays  a  sensitive  role  in  the  penetration  process.  In  other 
words,  most  of  the  energy  remains  in  the  projectile  as  it  passes  through  the 
target  only  a  small  portion  of  which  is  used  to  damage  the  material  surrounding 
the  hole.  If  this  difference  between  the  full  scale  and  sub-scale  system  is  to 
be  neglected,  then  scaling  is  no  longer  the  issue. 
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